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Abstract 

We explore the homogenization limit and rigorously derive upscaled equations for 
a microscopic reaction-diffusion system modeling sulfate corrosion in sewer pipes 
made of concrete. The system, defined in a periodically-perforated domain, is semi- 
linear, partially dissipative and weakly coupled via a non-linear ordinary differen- 
tial equation posed on the solid-water interface at the pore level. Firstly, we show 
the well-posedness of the microscopic model. We then apply homogenization tech- 
niques based on two-scale convergence for an uniformly periodic domain and derive 
upscaled equations together with explicit formulae for the effective diffusion coeffi- 
cients and reaction constants. We use a boundary unfolding method to pass to the 
homogenization limit in the non-linear ordinary differential equation. Finally, be- 
sides giving its strong formulation, we also prove that the upscaled two-scale model 
admits a unique solution. 

Key words: Sulfate corrosion of concrete, periodic homogenization, semi-linear 
partially dissipative system, two-scale convergence, periodic unfolding method, 
multiscale system. 



1 Introduction 

This paper treats the periodic homogenization of a semi-linear react ion- diffusion 
system coupled with a nonlinear differential equation arising in the modeling 
of the sulfuric acid attack in sewer pipes made of concrete. The concrete cor- 
rosion situation we are dealing with here strongly influences the durability of 
cement-based materials especially in hot environments leading to spalling of 
concrete and macroscopic fractures of sewer pipes. It is financially important 
to have a good estimate on the moment in time when such pipe systems need 
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to be replaced, for instance, at the level of a city like Los Angeles. To get 
good such practical estimates, one needs on one side easy-to-use macroscopic 
corrosion models to be used for a numerical forecast of corrosion, while on the 
other side one needs to ensure the reliability of the averaged models by allow- 
ing them to incorporate a certain amount of microstructure information. The 
relevant question is: How much of this oscillatory-type information is needed 
to get a sufficiently accurate description of the heterogeneous reality? Due to 
the complexity of possible shapes of the microstructure, averaging concrete 
materials is far more difficult than averaging metallic composites with rigor- 
ously defined well-packed structure. In this paper, we imagine our concrete 
piece to be made of a periodically-distributed microstructure. Based on this 
assumption, we provide here a rigorous justification of the formal asymptotic 
expansion performed by us (in [TJ) for this reaction-diffusion scenario. Note 
that in [I] upscaled models are derived for a more general situation involving 
a locally-periodic distribution of perforation^]. Locally periodic geometries 
refer to a special case of x-dependent microstructures, where, inherently, the 
outer normals to (microscopic) inner interfaces are dependent on both spatial 
slow variable, say x, and fast variable, say y. 

In the framework of this paper, we combine two-scale convergence concepts 
with the periodic unfolding of interfaces to pass to the homogenization limit 
(i.e. to £ 0, where e is a small parameter linked to the relative size of the 
perforation) for the uniformly periodic case. Here, the outer normals to the in- 
ner interfaces are dependent only on the spatial fast variable. For more details 
on the mathematical modeling of sulfate corrosion of concrete, we refer the 
reader to [2ll3] (a moving-boundary approach: numerics and formal matched 
asymptotics) , [I] (a two-scale react ion- diffusion system modeling sulfate cor- 
rosion), as well as to [5], where a nonlinear Henry-law type transmission con- 
dition (modeling H 2 S transfer across all air-water interfaces present in this 
sulfatation problem) is analyzed. Mathematical background on periodic ho- 
mogenization can be found in e.g., [SHE], while a few relevant (remotely re- 
sembling) worked-out examples of this averaging methodology are explained, 
for instance, in [9|10|I11|I12|I13|I14] . It is worth noting that, since it deals with 
the homogenization of a linear Henry-law setting, the paper [UJ is related to 
our approach. The major novelty here compared to [UJ is that we now need 
to pass to the limit in a non-dissipative object, namely a nonlinear ordinary 
differential equation (ode). The ode is describing sulfatation reaction at the 
inner water-solid interface - place where corrosion localizes. This aspect makes 
a rigorous averaging challenging. For instance, compactness-type methods do 
not work in the case when the nonlinear ode is posed on e-dependent surfaces. 
We circumvent this issue by "boundary unfolding" the ode. Thus we fix, as 
independent of e, the reaction interface similarly as in [15], and only then we 
pass to the limit. Alternatively, one could use varifolds (cf. e.g. [H]), since this 



The word "perforation" is seen here as a synonym for "pore" or "microstructure". 
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seems to be the natural framework for the rigorous passage to the limit when 
both the surface measure and the oscillating sequences depend on e. However, 
we find the boundary unfolding technique easier to adapt to our scenario than 
the varifolds. 

Note that here we approach the corrosion problem deterministically. However, 
we have reasons to expect that the uniform periodicity assumption can be 
relaxed by assuming instead a Birkhoff-type ergodicity of the microstructure 
shapes and positions, and hence, the natural averaging context seems to be 
the one offered by random fields; see ch. 1, sect. 6 in [17] . ch. 8 and 9 in [T8] , 
or [19]. But, methodologically, how big is the overlap between homogenizing 
deterministically locally-periodic distributions of microstructures compared to 
working in the random fields context? We will treat these and related aspects 
elsewhere. 

The paper is organized as follows: We start off in section [2] (and continue 
in section [3]) with the analysis of the microscopic model. In section HI we 
obtain the e-independent estimates needed for the passage to the limit e — > 0. 
Section [5] contains the main result of the paper: the set of the upscaled two-seal 
equations. 

2 The microscopic model 

In this section, we describe the geometry of our array of periodic microstruc- 
tures and briefly indicate the most aggressive chemical reaction mechanism 
typically active in sewer pipes. Finally, we list the set of microscopic equa- 
tions. 

2.1 Basic geometry 

Fig. [TJ (i) shows a cross-section of a sewer pipe hosting corrosion. We assume 
that the geometry of the porous medium in question consists of a system of 
pores periodically distributed inside the three-dimensional cube Q := [a, b] 3 
with a, b G R and b > a. The exterior boundary of Q consists of two disjoint, 
sufficiently smooth parts: T N - the Neumann boundary and T D - the Dirichlet 
boundary. The reference pore, say Y := [0, l] 3 , has three pairwise disjoint 
connected domains Y s , Y w and Y a with smooth boundaries T sw and T wa , as 
shown in Fig. [1] (iii). Moreover, Y := Y s U Y w U Y a . 

Let e be a sufficiently small scaling factor denoting the ratio between the 
characteristic length of the pore Y and the characteristic length of the domain 
Q. Let x w an d X a De the characteristic functions of the sets Y w and Y a , 
respectively. The shifted set Y™ is defined by 

Y™ :=Y + Z 3 j=0 k jej for k = (k u k 2 , k 3 ) G Z 3 , 

where ej is the j th unit vector. The union of all shifted subsets of Y™ multiplied 
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Fig. 1. Left: Cross-section of a sewer pipe pointing out one region. Middle: Periodic 
approximation of the periodic rectangular domain. Right: Reference pore configu- 
ration. 

by e (and confined within Q) defines the perforated domain Q £ , namely 

n £ := U fceZ s{ei^ | sY™ C Q}. 

Similarly, Q\, T S E W , and Y™ a denote the union of the shifted subsets (of fl) 
Y fe a , T™, and T^ a scaled by e. Since usually the concrete in sewer pipes is not 
completely dry, we decide to take into account a partially saturated porous ma- 
teriaO} We assume that every pore has three distinct non-overlapping parts: 
a solid part (grain) which is placed in the center of the pore, the water film 
which surrounds the solid part, and an air layer bounding the water film and 
filling the space of Y as shown in Fig. [TJ The air connects neighboring pores to 
one another. The geometry defined above satisfies the following assumptions: 

(1) Neither solid nor water-filled parts touch the boundary of the pore. 

(2) All internal (air-water and water-solid) interfaces are sufficiently smooth 
and do not touch each other. 

These geometrical restrictions imply that the pores are connected by air-filled 
parts only which is needed not only to give a meaning to functions defined 
across interfaces, but also to introduce the concept of extension as given, for 
instance, in [20J. Furthermore, there are no solid-air interfaces. 

2.2 Description of the chemistry 

There are many variants of severe attack to concrete in sewer pipes, we focus 
here on the most aggressive one - the sulfuric acid attack. The situation can 
be described briefly as follows: (The anaerobic bacteria in the flowing waste 
water release hydrogen sulfide gas (H 2 S) within the air space of the pipe. These 
bacteria are especially active in hot environments. From the air space inside 
the pipe, H 2 S(gjiE\ enters the pores of the concrete matrix where it diffuses 
and then dissolves in the pore water. The aerobic bacteria catalyze some of the 
H 2 S into sulfuric acid H 2 SO^. H 2 S molecules can move between air-filled part 
and water- filled part the water-air interfaces [21]. We model this microscopic 

2 The solid, water and air parts corresponds to Y s , Y w and Y a , respectively. 

3 H 2 S(g) and H 2 S(aq) refer to gaseous, and respectively, aqueous H 2 S. 
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interfacial transfer via Henry's law [22], (see the boundary conditions at T^ 0- 
in (J3J) and (@|). H 2 S0 4 being an aggressive acid reacts with the solid matrbo 
at the solid-water interface, which is made up of cement, sand, and aggregate, 
and produces gypsum (i.e. CaSO^ ■ 2H 2 0). Here we restrict our attention to 
a minimal set of chemical reactions mechanisms as suggested in [2] , namely. 

' 10H + + SOl 2 + org. matter — ► H 2 S(aq) + AH 2 + oxid. matter 
H 2 S{aq) + 20 2 — > 2H + + SOf 
H 2 S(aq) ^ H 2 S(g) 
2H 2 + H + + SOf + CaC0 3 — ► CaS0 4 ■ 2H 2 + HCO^ 

(1) 

We assume that reactions ([T]) do not interfere with the mechanics of the solid 
part of the pores. This is a rather strong assumption since it is known that (TTJ) 
can actually produce local ruptures of the solid matrix [23J . For more details on 
the involved cement chemistry and connections to acid corrosion, we refer the 
reader to [24J (for a nice enumeration of the involved physicochemical mecha- 
nisms), [23] (standard textbook on cement chemistry), as well as to [23|26II27] 
and references cited therein. For a mathematical approach of a similar theme 
related to the conservation and restoration of historical monuments, we refer 
to the work by R. Natalini and co-workers (cf. e.g. [28]). 



2. 3 Setting of the equations 

The data and unknown are given by 

u\ Q : Q — > R + - initial concentration of H 2 SOi{aq) 
u 20 : Q — > R + - initial concentration of H 2 S(aq) 
u £ 30 : Q — > W + - initial concentration of H 2 S(g) 
u\ : Q — > R + - initial concentration of moisture 
«5 : Q — > R + - initial concentration of gypsum 

tif : Tjj X (0,T) — > R + - exterior concentration (Dirichlet data) of H 2 S(g) 
ul:tt £ x (0,T) — > R - concentration of H 2 S0 4 (aq) 
u 2 :Q £ x (0,T) — y R - concentration of H 2 S(aq) 
u £ 3 -.nix (0, T) — > R - concentration of H 2 S(g) 
u\ : Q £ x (0, T) — > R - concentration of moisture 
u\ : Tf" x (0,T) — )• R - concentration of gypsum 

All concentrations are viewed as mass concentrations. We consider the fol- 
lowing system of mass-balance equations defined at the pore level. The mass- 



4 The solid matrix is assumed here to consist of CaCO% only. This assumption can 
be removed in the favor of a more complex cement chemistry. 
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balance equation for H2SO4 is 

d t u\ + div{-d\Vu\) = -k\u\ + k £ 2 u £ 2 , x G Q £ , t G (0, T) 
Ui(x,0) = u £ 10 (x), xeQ £ 
-n £ ■ d\Wu\ = 0, x G rf\ £ G (0, T) U 
-n £ ■ d\Vu\ = er]{ul, u%), x G q w £ G (0, T). 

The mass-balance equation for H 2 S(aq) is given by 

d t u £ 2 + div(-d £ 2 Vu £ 2 ) = k\u\ - k £ 2 u £ 2 , x G f2 e , £ G (0, T), 
«IO,0) = 4 (x), x G fi £ 
-n e • d £ 2 Vu £ 2 = e{a £ {x)ul - b £ (x)u £ 2 ), x G Tf a , t G (0, T) U 
-n e • d £ 2 Vu £ 2 = 0, x G rf , £ G (0, T). 

The mass-balance equation for H2S(g) reads 

d t u £ 3 + div{-d £ 3 Vul) = 0, iGflj, £ G (0, T) 
4(x,0) = w| (x), i6fl; 
-n £ • c/|V^ = 0, x G T^, £ G (0, T) (4) 

4(x, £) = uf(x, £), x G T D , t G (0, T) 
-n £ ■ df,Vu e 3 = -e(a £ (x)u £ 3 - b £ (x)u £ 2 ), x G T™ a , £ G (0, T). 

The mass-balance equation for moisture follows 

d t u\ + div(-dlVul) = k{u\, x G tt £ , £ G (0, T) 
t4(x, 0) = «4 (x), 

-n e • d*V< = 0, x G r™, £ G (0, T) U 
-n £ ■ d\Vu\ = 0, x G rf , £ G (0, T). 

The mass-balance equation for the gypsum produced at the water-solid inter- 
face is 

d t ul = V (ul,ul), XGIT, te(0,T) 
u £ 5 (x,0) = ut (x), x G Tf", £G(0,T). 1 } 



3 Weak formulation and basic results 

We begin this section with a list of notations and function spaces. Then we 
indicate our working assumptions and give the weak formulation of the mi- 
croscopic problem; we bring reader's attention to the well-posedness of the 
system 

3. 1 Notations and function spaces 

We use (a,/3)( ,T)xfF := Jo /n« a/3dxdt, (a, /3) (0 ,T)xr e := Jo Jr e a(3da x dt. (•), 
I • I and || • || denote the dual pairing of iJ 1 (fi e ) and H -1 ^ 6 ), the norm in 
L 2 (Q £ ), and the norm in f/" 1 (f2 e ), respectively. tp + and tp~ will point out the 
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positive and respectively the negative part of the function ip. We denote by 
C#{Y), H#(Y), and H^(Y)/M., the space of infinitely different iable functions 
in MJ 1 that are periodic of period Y, the completion of CT(Y) with respect 
to if 1 — norm, and the respective quotient space, respectively. Furthermore, 
H* D (Sl) := {u G H\9)\u = on T D }. The Sobolev space H P (Q) as a com- 
pletion of C°°(Q) is a Hilbert space equipped with a norm 

and (cf. Theorem 7.57 in the embedding H 13 ^) ^ L 2 (f2) is continuous. 
Since we deal with an evolution problem, we need typical Bochner spaces like 
L 2 (0, T; if 1 (!])), L 2 (0, T; L 2 (0, T; #^(0)), and L 2 ((0, T) xft; 

In the analysis of the microscopic model, we use frequently the following trace 
inequality for e— dependent hypersurfaces T™ a : For ip E G H 1 ^ 8 ), there exists 
a constant C*, which is independent of e, such that 

e\(Pe\h { r s) < C*(\fs\h m + £ 2 \V<fe\h m )- (7) 

The proof of ((Tj) is given in Lemma 3 of [30]. For a function y? e G H 13 ^ 6 ) with 
/3 G (5, 1), the inequality ((7j) refines into 

1 12 12 . 2/3 r r i^ £ ( x ) -v e {y)\ 2 i , x ,m 

eb e |x»(r.) < c (|^L w + £ [a: _ y | W+2/3 ( 8 ) 

where Cq is again a constant independent of e. For proof of (jHJ), see [15] . 
To simplify the writing of some of the estimates, we employ the next set of 
notations: 



£ I 

i I ' 



di := min_ I d e f I, i 6 {1,2,3,4}, di := min_ I d, 

[0,T]xn [0,T]xQ 

D m := max l&dfL I, m G {1, 2, 3}, := min I /cf I, j G {1,2} 

[o,nxn« m J [o,T]xn L J 

X,- := min_ I I, kj := min_ I /cf I, 

[o,T]xn J [o,r]xn J 

/c m . = sup I /c m I , /c m . = sup I k m I , 

(0,T)xH (0,T)xQ 

K%:= sup l^/c^l, M<:= sup | u\ |, z G {1, 2, 3, 4, 5}, 

(0,T)xfi (0,T)xn 





:= sup |a e |, B°° :- 
(o,r)xr™ a 


= sup |6 e |, 
(o,T)xr™ a 




:= sup |<9 t a e |, B°° 
(o,T)xr™ a 


:= sup |9 t fe e |, 
(o,T)xr™ a 




:= sup a , 6°° := 

(o,T)xr u,a 


sup |6|, 
(o,T)xr™ a 


g°° 


:= sup \Q(s)\, 

sG(0,T)xr|™ 


v-= WvWoo, v-= 



\\dtV\ 
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3.2 Assumptions on the data and parameters 



We consider the following restriction on the data and parameters: 

(Ai) d t e l°°((o,t) x y) 3x3 , g l°°((o,t) x y) 3x3 , d tt di e l°°((o,t) x 

F) 3x3 , (di(t,x)&£) > d i0 | C | 2 for d i0 > 0, for every £ G M 3 , (t.z) G 
(0,T)xy,i 6 {1,2,3,4}. 
(A2) 77 is measurable w.r.t. t and x and r/(a, /3) = fc|-R(o;)Q(/3), is sub-linear 
and locally Lipschitz function and Q is bounded and locally Lipschitz 
function such that 

positive, if a > 0, positive, if /3 < (3 max , 

R(a) = i Q{(3) = i 

( 0, otherwise ( 0, otherwise 

Additionally to (A2), we sometimes assume (A2)', that is 
(A2)' d t r] < fj. 

(A3) ufo G L 2 (ft £ ) n i G {1,2,4}, u% G L\n\) fl L-(Of), u§ G 

L 2 (rf ) ni£(rf )■ 

(A4) a°°M 3 = 6°°M 2 , fc^Mi = M 4 , fcjMi = fcfM 2 . 

(A5) a,be C 1 ([0,T];C ' a! (r wa )),a,6> in [0,T] x T wa , d t a, d t b G L°°((0,T) x 
(A6) c^uf, d tt uf and Vd t u® are bounded. 

(A7) A; 3 G C rl ([0,^;C°' a (r aw )) and kj G C 1 ([0, T]; C ' a (F)) for any j G {1,2} 
and a g]0, 1]. 

The assumptions (Al)-(A3), (A5), and (A6) are of technical nature. The first 
equality in (A4) points out an infinitely fast (equilibrium) Henry law, while the 
last two equalities remotely resemble a detailed balance in two of the involved 
chemical reactions. 

3. 3 Weak formulation of the microscopic model 

Definition 1 Assume (Al) and (A3). We call the vector u £ = («f , u|, u|, w|, 
a weak solution to Qty-fflifu] G L 2 (0,T; H l {Q, e )),d t u e G L 2 (0,T; i?~ 1 (f2 e )),j G 
{1,2,4}, u% G u° + L 2 (Q,T;Hi D ({ll)),d t ul G uf + L 2 (0, T; ff-^flf)),^ G 
L°°((0,T) x rf°),5^«| G L°°((0,T) x rf ) such that the following identities 
hold 

(d t u\, <pi)(o,T)xa° + (diVuf), Vy?i)( ,T)xfF 

= -(^1^1, </ ? i)(o,T)xn« + (^2^2, </ ? i)(o,T)xn« (9) 
- e(»7(«i J «H) J ¥'i)(o,r)xrf» J 
(d t ul, tp 2 ) (o,T)xne + (d|Vw|), V</? 2 )(o,T)xn« 

= (*i u l> ^2)(0,r)x^ - (^2 M 2 ? ¥>2)(0,T)xn« (10) 

+ e(a e itf, <£>2)(o,T)xr™» — s(a £ u £ 2 , v ? 2)(o,T)xr™«, 
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(d t u £ 3 , ^ 3 }(o,T)xnf = -(d £ 3 Vu £ 3 ), V</? 3 )(o,T)xnf ^ 
— e(a £ M3, V3)(o,T)xr™ a + ^(o, £ u £ 2 , v?3)(o,T)xr™a, 

(^tt4> V ? 4)(0,T)x^ = -(t^V«|), V<p 4 )(0,T)xn* + V ? 4)(0,T)x^ (12) 

for all ipj E L 2 (0,T;H\tt £ )),j E {1,2,4} and ip 3 E L 2 (0, T; flJ D (fif)) to- 
gether with the ode 

d t u\= a.e. on (0,T) xTf s (13) 

and i/ie initial conditions 

u £ (0, x) = u £ (x) x Eft 6 for all i E {1,2, 4}, 
u £ 3 (0,x)=ul {x) xEQl, (14) 



5.^ Basic results 

Lemma 2 (Positivity and L°° -estimates) Assume (A1)-(A6), and let t E 
[0, T] 6e arbitrarily chosen. Then the following estimates hold: 

(i) u £ (t) > 0, « G {1,2, 4} a. e. in Q £ , u £ 3 (t) > a. e. Q\ and u e 5 (t) > a. e. on 

e ■ 

(ii) u £ {t) <Mi,iE {1,2}, u\{t) < (t+l)M 4 a.e. mVt £ , u £ 3 {t) < M 3 a.e. in Of 
and u|(t) < M 5 a.e. on r™ s . 

Proof (i). We test fl9l- ffl~2|) with cp = (— nf _ , — u 2 ~ , — u £ 3 ~ , — «|~) element of 
the space [L 2 (0, T; tf 1 ^'))] 2 x ^ 2 (0, T; ^(fif) x L 2 (0, T; H l {Q, e ). We obtain 
the following inequality 

l -d t \u £ ~\ 2 + rfilVnM 2 < -k x \u\-\ 2 + fc^K", (15) 

-e{r){ul,ul),-u\-) rfU . 



Note that the first term on the r.h.s of ffl5|) is negative, while the third term 
is zero because of (A2). We then get 

d t \u\-\ 2 + 2di\Vu\~\ 2 < [\u\-\ 2 + |nr| 2 ) • (16) 

On the other hand, (fTUI) leads to 



\d t \ u r\ 2 + d 2 \vu £ 2 -\ 2 < h \ (k-| 2 + \u\-\ 2 

+ £a°°(u £ 2 ~ ,u £ 3 ~) r ^ + eb°°\u £ 2 '\lu,a. 
By the trace inequality ([7j) (with e < 1), we get 



d t |n|~| 2 + 2(rf 2 - C*6°°)|Vnr| 2 < fc?° (K~| 2 + 
+ 2C*o°°K~| 2 + 2ea°°(u £ 2 ',u £ 3 ') rr . 



ur\ 2 



(17) 
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(TTTj) leads to 

d t \u £ f\ 2 + 2{d 3 - C*a°°)\Vu £ f\ 2 < 2eb co {u £ 2 -,u £ 3 -) Tr + 2CV>r| 2 , (18) 
while from (fT2|) . we see that 

d t \ui~\ 2 + 2d 4 |v«ri 2 < &r (k~i 2 + wr\ 2 ) . (19) 

Adding up inequalities ( TIE]) - (TI51) gives 



i=i 



and hence, 



otJICI 2 + 2di|Vur| 2 + 2(d 2 - C*b°°)\Vu £ 2 -\ 2 

+ 2(d 3 - C*a°°)\Vu £ 3 -\ 2 + 2rf 4 |Vur| 2 

4 

< (2A^° + + 2C*b°° + 2C*a°°) £ 
+ 2e(a~ + 6~)(ii|~ itf-)r»», 



^ E K'l 2 + 2rfi|V«ri 2 + 2 (rf 2 - C*b°°)\Vu £ 2 -\ 2 



(20) 



1=1 



+ 2(d 3 - C*a°°)|Vur| 2 + 24|Vm £ -' 2 



< (2k™ + k™ + C*(a°° + 6°°)) E 



i,' - ( 21 ) 



l<-| 2 



8=1 



+ e (a°° + &°°)(%rir- + ^"Ir™ 

Applying the trace inequality ([7]) to estimate the last term on the right side 
of (I2T]) . we finally get 



$ E K"! 2 + 2rfi|V«r| 2 + (24 - 2C*6°° - C*5{a°° + b°°))\Vu 

i=l 

+ (2d 3 - 2C*a°° - ^(a°° + &°°))|Vwr| 2 + 2d 4 |V^ 


4 



e-|2 
2 I 



<^E 



\uT?. 



1=1 

Thus, we have 

i=l i=l 

where C x := 2k^ + k^ + C*(a°° + 6°°) + C*(5 + ±)(a°° + 6°°) and 5 is chosen 
conveniently Gronwall's inequality together with [wf(0)] _ = gives now the 
desired result. Note that (A2) ensures automatically the positivity of u%. 
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(ii). We consider the test function 

(<p 1} <p 2 , <p 3 , <f 4 ) = (K - K - M 2 )+ K - M 3 ) + , « - {t + 1)M 4 )+). 

Obviously, G [^(OJjiJ 1 ^))] 2 x L 2 (0, T; i^ D (fif) x L 2 (0,T; H 1 ^) is 
allowed as test function. We obtain from (jUJ) that 

^|K - Mx) + | 2 + d!|V(tt? - Mi) + | 2 < -h\(ul - M0 + | 2 

- (fciMj, K - M0+) 

+ A^((<-M 1 )+(u§-M 2 ) + ) 

+ (A; 2 °°M 2 , K - M0+) 

-£(r ? ( M f, M £ 5 ),K-M 1 ) + ) rr . 

Relying on (A4), we get the estimate 

d t \{u\-M x )+\ 2 < fc 2 00 (|K-M 1 )+| 2 + |(w|-M 2 )+| 2 ). (22) 
( ITU]) in combination with (A4) gives that 



d t \{u\ - M 2 )+| 2 + 2{d 2 - C*b°°)\V{u £ 2 - M 2 )+| 2 

<*r(i(«i-^i) + i 2 + i(«i-^) + i 2 j 

+ 2C*6°°|(m|-M 2 ) 



\2 



WO . 



(23) 



+ 2ea°°((u £ 2 - M 2 )+ K - M 3 )+) r 

By ( ITT]) , we obtain 

d t \(ul - M 3 ) + | 2 + 2(4 - C*0|V(4 - M 3 ) + | 2 

< 2C*a°°|V(^-M 3 ) + | 2 (24) 
+ 2 £ &°°(K - M 2 ) + , (ttf - M 3 )+) rr • 

Using again (A4), ( fT2]) yields 

- (* + 1)M 4 )+| 2 < k?(\(ut - M0 + | 2 + |K - (t + 1)M 4 )+| 2 ). (25) 
Adding up (l2"2"]) - (l25]) side by side, we get 

jZdt^u] - M,) + | 2 + d t \(ul - (t + 1)M 4 ) + | 2 + (2d 2 - 2C*b°°)\V(u £ 2 - M 2 ) + | 2 
+ (2rf 3 - 2C*a°°)|V(4 - M s) + | 2 

3 

< (2A£° + k™ + 2C*a°° + 2C*b°°)(Y / ~ Mj) + \ 2 

3=1 



+ IK - (t + 1)M 4 )+| 2 ) + e(a°° + b°°)(8\(ul - M 2 )~ 
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+ ^|K-M 3 ) + | 2 r )- 
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We use the trace inequality ([7]) (with e < 1) to deal with the boundary terms in 

(1261) . Then Gronwall's inequality yields for all t G (0, T) the following estimate 
u){t) < Mj, j G {1, 2, 5} a. e. in Sl s 

u £ 3 (t) < M 3 , a. e. in Q[ 

u\<{t + 1)M 4 a.e. in Vt £ . 

Furthermore, by (A2) u\ is bounded. 



Proposition 3 (Uniqueness) Assume (Al)-(A^). Then there exists at most 
one weak solution in the sense of Definition^ 



Proof. We assume that v?' 6 = (wi' £ , u 2 £ , u J 3 ' e , u^ £ , u J 5 ' E ), j G {1,2} are two dis- 
tinct weak solutions in the sense of Definition [TJ We set u\ := u\ £ — u 2,£ for 
all i G {1,2,3,4}. Firstly, we deal with (15). We obtain 

d t ul' £ - d t u 2 5 ' £ = r](u\' £ , u\' £ ) - n(ul' £ , u 2 5 ' £ ). (26) 

Integrating ( 126]) along (0,T) and using (A2), we get 

\u 5 ' £ - u/\ < kfc R c Q M l f \uY - ul> £ \dT + k?c R Q°° f \u{' £ - ul' £ \dr. 

Jo Jo 

Gronwall's inequality implies 

\u l /(t) - u/(t) \ < C 2 f \u{ £ - u{ £ \dr for a.e. t G (0,T), (27) 

Jo 

where C 2 := k^c R Q°°(l + C 3 te Cst ) and C 3 := k\ x> c R c Q M 1 . We calculate 

^d t \u\\ 2 + d 1 \Vu £ 1 \ 2 < -hint] 2 + k?{u\,u £ 2 ) + - r/ 2 ,uf)r r , (28) 

where we denote rji — n 2 := n(u\ ,£ , u\' £ ) — i](ul' £ , u^ £ ). We can write 

1 k°° 

-d t \u £ \ 2 + d^Vutf < -hluW 2 + ^-(K| 2 + K| 2 ) 

+ eC,{uY-u¥,u\)r r (29) 
+ ek™c R Q co (u{' £ -u 2 f,u\) Tf ». 



Now, inserting (12?]) in fl29|) yields 



I LOO 

-d t \u £ \ 2 + d^Vull 2 < -hluH 2 + -J-(K| 2 + \u £ 2 \ 2 ) 

eCl ft 



s~i I £\1 ^ 3 f i ei2 7 

~|~ Ls^E Mi \psw ~\~ - / Mi r™(IT. 

25 Jo E 



(30) 
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where C4 := kfcnQ™ + . Using (jTJ), we estimate the last two terms in f l30l) 
to obtain the inequality 

1 k°° 

-d t \u\\ 2 + rfiiv^i 2 < -hiuii 2 + -^(Kl 2 + \u £ 2 \ 2 ) + C*C A (\u\\ 2 

C 2 * ^ 



+ e 2 \Vul\) + C*^ J Q {\u\\ 2 + s 2 \Vu\\ 2 )dr. 



Note that the constant C*, arising from in (13T|) . stems from ([7j). Rearranging 
now the terms, we have 

d t \ul\ 2 + (2d, - 2C*C 4 e 2 )\Vu £ 1 \ 2 + 2k 1 \u £ 1 \ 2 < (A£° + C*C 4 )(K| 2 

C 2 ^ (32) 



/ {\u\\ 2 + e 2 \Vu\\ 2 )dr. 
Jo 

Following the same line of arguments as before, we obtain from fllUp that 

(33) 



d t \u £ 2 \ 2 + 2d 2 \Vu £ 2 \ 2 < -2k 2 \u e 2 \ 2 + k?(\u\\ 2 + u e 2 \ 2 ) 
+ 2£a 00 ( y u £ 3 ,u £ 2 ) rr + 2eb°°\u £ 2 \ 2 r , 



while from ffTTj) . we deduce 

d t \u £ 3 \ 2 + 2d 3 \Vu £ 3 \ 2 < 2eb QO (u £ 2 ,u £ 3 ) r ^ +2ea°°\ul\l r . (34) 

Proceeding similarly, (fT2"|) yields 

d t \ul\ 2 + 2rf 4 |V<| 2 < A; 2 °°(K| 2 + \u\\ 2 ). (35) 
Putting together (I3"2"|) -fl3"5l). we get 



d t Y>t=M\ 2 + (2di - C*C 4 £ 2 )|V^| 2 + 2d 2 |Vu|| 2 + 2rf 3 |V^ 12 

E|2 , nl. I„.£|2 



+ 2rf 4 |Vw|r + 2fci|wf 1 
< (2k™ + k™ + C*C 2 )Zt =1 \u £ \ 2 

+ C*^ J\\u £ \ 2 + e 2 \Vu £ \ 2 )dT 
+ 2eb\u 2 \ 

+ e(a°° + b°°)(5\u £ 2 \ 2 rwa + ~|uf |£»„). 

Applying the trace inequality ([7]) to the boundary terms in f[3"6"j) . we get 

a t Sti|<| 2 + (2d! - 2C*C 4 e 2 )\Vul\ 2 

+ (2d 2 - 2C*b°°e 2 - C*fe 2 (a°° + b°°))\\7u £ 2 \ 2 

+ (2d 3 - 2C*a°°e 2 - ^^(a°° + 6°°))|Vm|| 2 



+ 2d A \Vu\\ 2 + 2A; 1 |^| 2 < C 5 £i 1 KI 2 



(36) 
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+ C*^ j\\u\\ 2 + e 2 \Vu{\ 2 )dr, 



(37) 



where C 5 := 2k™ + kf + C*C 2 + 2C*(a°° + b°°) + C*(a°° + b°°)(5 + ]). Let us 
choose £ and 5 such that 



e E 
5e 



2d! 



C*e 2 (a°° + b°°) 2d 2 - C*b°°e 2 



2d 3 - C*a°°e 2 ' C*e 2 (a°° + b°°^ 



With this choice of (e, 5), (!37|) takes the form 



d t Zt =1 \ul\ 2 + C\Vul\ 2 + C\ul\ 2 <C 6 (Zt=M\ 2 + / (\u\\ 2 + e 2 \Wu\\ 2 )dr) 

Jo 



where C 6 := 2A;~ + k™ + C*C 2 + C*(a°° + 6°°) + and C := min{2c/ 1 - 

2C*C26 2 , 2ki}. Taking in ( 137|) the supremum along t e (0,T) and applying 
Gronwall's inequality, we obtain the following estimate 



StiK | 2 + C / T |V^| 2 rft + C / T K| 2 dt < 0. 

JO JO 



(38) 



Thus, the proof of Proposition [3] is completed. 



Theorem 4 (Global Existence) Assume (Al) — (A3). Then there exists at 
least a global-in-time weak solution in the sense of Definition^ 



Proof. The proof is based on the Galerkin argument. Since the proof is rather 
standard, and here we wish to focus on the passage to the limit e — > 0, we 
omit it. 



4 A priori estimates for microscopic solutions 

This section includes the e— independent estimates. 

Lemma 5 Assume (A1)-(A6). Then the weak solution of the microscopic 
model (Ep-JT^j) satisfies the following a priori bounds: 



u) ||^(o,r^i(n«))< C, j E {1,2,3,4} (39) 

Vd t u £ 2 ||l2(o,t ; l2(^))< C, (40) 

d t Uj ||l2(0,T;L2(^))< C, (41) 

u 3 ||i 2 (0,T;J^ 1 (r2f))< C, (42) 

Vd t U £ 3 ||L2(0,T;L2(nf))< C, (43) 

d t u% \\mo,T;L\n\))< C, (44) 
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II u l IU°°((o,T)xri w )^ C) 
II d t ul ||L 2 ((o,T)xr|«')< C. 

In ^3^)-^J^, the generic constant C is independent of e. 

Proof. We test ([9]) with ipi = u\ to get 

-dt\u\\ 2 + d 1 \Vu\\ 2 < -/ci|^| 2 + fc 2 °°K,^) -e(r],ul)r™, 

k°° 

< -f(K\ 2 + \u e 2 \ 2 )+Ek™Q™ CR (ulu\)T r . 

After applying the trace inequality to the last term on r.h.s of (JSJ), we 

1 k°° 

-d t \u\\ 2 + d^VuW 2 < -J-(K| 2 + \u £ 2 \ 2 ) + C*k?Q°°c R (\ul\ 2 + e 2 \Vul\ 2 ) 

^d t \u\\ 2 + {d x - e 2 C*k™Q°°c R )\Vu\\ 2 < C 7 (\u\\ 2 + \u £ 2 \ 2 ), 
where C 7 := ^- + C*A;f Q°° Ci? . Taking y 9 = ^ in (flO]), we get 
1 

~<9 t |u|| 2 + d 2 \Vu £ 2 \ 2 < -^-{\u\\ 2 + |«§| a ) - k 2 \u £ 2 \ 2 
+ ea°°(u £ 3} u £ 2 )r^ +£&°°|«||r^a. 

Application of the trace inequality ([7]) only to the last term leads to 

1 k°° 

-d t \ul\ 2 + (d 2 - C*b°°e 2 )\Vu £ 2 \ 2 < ^(K| 2 + \ul\ 2 ) + 2 

We choose <p 3 = m| as a test function in (TTTT) to calculate 

^t|m 2 + (d 3 -C*a°V)|V?4| 2 < e6°°(u|,M|)r r + C*a°°m| 2 . 

Setting y9 4 = u\ in ( IT2|) . we are led to 

1 k°° 

^d t \ul\ 2 + d A \Vul\ 2 <^-{\u\\ 2 + \ul\ 2 ). 

Putting together (HHl) - (l5Tl) . we obtain 

istx^KI 2 + (dt - e 2 C*k™Q™c R )\Vu\\ 2 + d 4 |Vu 4 | 2 

+ (d 2 - C*6°V)|V^| 2 + {d 3 - C*a°°e 2 )|V^| 2 

k°° 

<(k? + ^- + C*b°° + C*a°°)^ =1 \u £ \ 2 
+ e(a°° + b°°)(u £ 3 ,u £ 2 ) rr . 
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Combing Young's inequality and the trace inequality to the boundary term, 
(152"]) turns out to be 

^t =1 d t \ut\ 2 + (di - e 2 C*k?Q°°c R )\Vul\ 2 

+ (d 2 - C*be 2 - ^-^(a°° + b°°))\Vu £ 2 \ 2 

+ (d 3 - C*ae 2 - ^Ha°° + b°°))\Vul\ 2 + d 4 |V<| 2 
2d 

LOO 1 

< (k? + -f + C*(a°° + &«)(* + ^M=iK\ 2 . 

Choosing e small enough and 5 conveniently such that the coefficients of the 
terms involving |Vwf | 2 and |Vw|| 2 stay positive, we are led to 



£tid t |wf| 2 + d'^VuW 2 + d' 2 \Vu £ 2 \ 2 + d' 3 \Vu £ 3 \ 2 + 2d 4 \Vu £ 4 \ 2 < C 7 E* =1 |u: 
where 

d[ := 2{d x - e 2 C*k™Q ( 

C*e 2 5 



d[ 


:= 2(4 


- E 2 C*k™ 




:= 2{d 2 


_ C * b oo £ 2 


4 


■ = 2(d 3 


- C*a°°e 2 



(a^ + fe 00 )), 



2 

C*hp 2 

^ r (« oo + n), 

while the constant C is given by 

LOO 1 

C 8 := 2k? + + C*a°° + + C*(a°° + 6°°)(5 + -). 

2 o 

Summarizing, we have 

I 2 + rfoS^lV^I 2 + d |Vm 2 < CEtiKI 2 , (53) 
where do := minfd'^ 4, gKj, 4}- By Gronwall's inequality, we have 



n=M\ 2 < cn =l \um\\ 

and hence 



II u £ j \\v(p,T;V(n>))< C for all i e {1,2,4} and ||w|IU 2 (o,T;£ 2 (^)) < C, (54) 

where C depends on initial data and model parameters but is independent of 
e. Integrating (j55|) along (0, T), we get 

II ||L2(0,T;Hi(^)) <C, j e {1,2,4}, 

II e II ^ /-« v* 3 ^ 

II M 3 ||£ 2 (o,T;H 1 (fif)) S 
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With the help of (A2) together with the boundedness of u\ , we conclude from 
([MD that 

II u l IU°°((o,r)xri m )< C. 
Multiplying ffl3|) by d t u% and using (A2), we get 



L 2 ((0,T)xr|™) 



< C. 



Now, we focus on obtaining e— independent estimates on the time derivative 
of the concentrations. Firstly, we choose <pi = d t u\ and get 

nd t u\d t u\dxdr + [ [ dTVuV\7 'd t u\dxdr 
j e Jo Jci s 

= — / k\u\d t u\dxdT + k £ 2 u £ 2 d t u\dxdT (56) 

Jo Jn s Jo Jn* 

— £ I / rjdtulda x dr. 
Jo Jv% w 

Consequently, it holds 

f f \d t u\\ 2 dxdr + f f (\d t {d\\Vu\\ 2 )-{d t d\)\Vu\f]dxdr 
Jo Jn E Jo Jn? \2 J 

< -7T f / d t \u\\ 2 dxdT 
2 Jo Jq. e 

(d t (r]ul) - (d t rj)u\)d<T x dT, 



rt 

io Jrz 



(1 _ *f^) /"* /" |^f| 2 fcdr <D 1 f [ \S7u\\ 2 dxdr 
2 jo Jn e Jo Jn? 

d°° { k°° /"* /" 

+ ^-/ |Vui | 2 dx + -pr- / / |?4| 2 Gfedr 
2 Jn e 2o Jo JfF 

+ |/ (\V\ 2 + K\ 2 + \V(0)\ 2 + K(0)\ 2 ) da x 

(57) 

where 77(0) := r/(ttf(0), w|(0)). Applying (J7J) and recalling ( )55l) . we have 



/* / |<9 t ^| 2 dxdr<C 9 , (5* 

JO JfK 



where 



C 9 :=I>i f / \Wu\\ 2 dxdT + — I \u\{Q)\ 2 dx+ ( ^- [ \Vu 10 \ 2 dx 
Jo J£F 2 JQ e 2 Jn? 
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+ 177 \ut\ 2 dxdr+ £ - [ (\fj\ 2 + \V(0)\ 2 + \€) 
2o Jo iff 2 jy% w 

+ ? f I (\u £ 1 \ 2 + s 2 \Vu £ 1 \ 2 )dxdr + ^- [ {\u\\ 2 + e 2 \Vu\\ 2 )dxdr, 
2 Jo iff 2 iff 



and 5 e 





w 5 £.00 
/v 2 



Testing ffTUj) with <p 2 = d t u 2 gives 



/ / \d t u £ 2 \ 2 dxdr + f I {-d t {d £ 2 \Vu £ 2 \ 2 ) - {d t d £ 2 )\Vu £ 2 \ 2 )dxdr 
Jo iff Jo iff 2 



K rt r h°° rt r 1 



<~T d t \u^dxdr+^- / / (^K 

Z JO iff 2 Jo iff 



2 

and hence, 

do 



+ 5|<9 t u e 2 | 2 )cfedT + I [ (|m 2 + \d t u £ 2 \ 2 ) da x dr 
2 Jo ir™ a v 7 

/ / <9 t |u 2 | 2 dcr-cdr, 
Jo Jrf a 



eb°° ft 



[ [ \d t u%\ 2 dxr+— [ \Vu 2 £ \ 2 dxdr 
Jo iff 2 iff 

\Vu £ 2 {0)\ 2 dx + D 2 f l f \Vu £ 2 \ 2 dxdr 

I s Jo iff 

f / (i|^| 2 + 5|^ £ 2 | 2 )rfxc/r 
Jo iff o 



< 2 



2 iff 



+ 2 
+ 



/ / (\u £ 3 \ 2 + e 2 \Vu £ 3 \ 2 + e 2 \Vd t u £ 2 \ 2 ) dxdr 



2 ir 

By ([7D and ([SSD, we get 



2 

+ ^rl (\u £ 2 \ 2 - \u £ 2 (0)\ 2 )da x . 



' C*a°° k?5^ 



ft \d t u £ 2 \ 2 dxdr<C w (l + e 2 f [ \V d t u £ 2 \ 2 dxdr 
Jo iff V Jo iff 



Consequently, choosing 6 e]0, 2 ^L°°° [, we are led to 



/* / \d t u £ 2 \ 2 dxdr < C 10 (l + e 2 f f \Vd t u £ 2 \ 2 dxdr), (59) 
Jo iff Jo iff 



where 



C W :=D 2 / |V^| 2 cfec2T + ^ 2 - / \Vu £ 2 (0)\ 2 dx + ^- / |^| 2 rfxdr 
io iff 2 iff 2o io iff 

+ J ne (\u £ 2 \ 2 + e 2 \Vu £ 2 \ 2 + \u £ M\ 2 + e 2 \Vu £ M\ 2 ) dx 

C*a°° r* 



+ 



f f {\ul\ 2 + e 2 \Vul\ 2 )dxdr. 

Jo iff 
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The initial data -u| holding in Q\ and the Dirichlet data u 3 acting on the 
exterior boundary of Vt\ are considered here as restrictions of the respective 
functions defined on whole of Q. Testing now ( II ip with ip 3 = dt(u 3 — u 3 ) leads 
to 



f ( \d t ul\ 2 dxdr + % ( \Vu 
Jo Jnc 2 Jn^ 



£ |2 

3 1 



< d il \Vu%m 2 + \(\dtu%\* + \d t u^) 
2 Jn^ 2 

+ D 3 f f |W 3 | 2 + ^/7 (\Vu £ 3 \ 2 + \Vd t uf 
Jo Jn? 2 Jo Jn e 

+ ^/7 K 3 r + y(«°° + ff \dtul\ 
o Jo Jrf a 2 Jo Jvf a 

2 Jo Jv™ a o Jo Jt 



u%\ 2 . 



Using (J7J) and (A6), we obtain 

/ / \d t u £ 3 \ 2 dxdr <C 11 (l + e 2 5 [ \Vd t u £ 3 \ 2 dxdr), (60) 
Jo Ju s Jo Jn s 

where 5 e]0, c , (a i +boo) [ and 

Cn := D 3 f f \Vu £ 3 \ 2 dxdr + ^ f |W 3 (0)| 2 cfe + 1 ff \Vu°\ 2 
Jo Jn~ 2 Jn? 2d Jo Jn e 

d°° ft r CI*/! 00 rt r 

+ %/ / (|V^ 2 + |V^f| 2 ) + ^/ / (\ul\ 2 + s 2 |W 3 
2 Jo Jn e o Jo Jn e 



+ / {\u £ 2 \ 2 + e 2 \Vu £ 2 \ 2 )dxdr 



5 JO JQ' 

+ C * {a °° n + b00) f f (\d t uf\ 2 + e 2 \Vd t uf\ 2 )dxdr. 



2 jo Jn 
From (TT21 . we get 



ff \d t u\\ 2 dxdT<C 12 . (61) 

Jo Jci e 

In order to estimate (1591) and (160]) . we proceed first with differentiating (flQl) 
with respect to time and then testing the result with d t u 2 . Consequently, we 
derive 



-/ \d t u £ 2 \ 2 dx + d 2 f I \Vd t u 2 £ \ 2 dxdr (62) 
2 Jn^ Jo Jn? 
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+ 



f I (\{d t d 2 \Vu 2 £ \ 2 -{d t d t d 2 )\Vu 2 £ \ 2 
Jo jj> \2 



<^/7 (i^i 2 +i^i 2 ) + ^/7 (ki 2 +i^i 2 ) 

z Jo jq e z Jo Jn e 



4I 2 



-fa I I \d t u £ 2 \ 2 - — ( ( d t \u 
Jo Jn e 2 Jo Jn s 

Ffl°° ft f 1 

+ (-\d t u £ 2 \ 2 + 5\d t u £ 3 \ 2 )da x dr 

2 Jo JVf a 

+ / \d t u £ 2 \ 2 da x dr + — / / d t \u £ 2 \ 2 da x dr. 

2 JO JT™ a 2 Jo JTf a 

Using ([7]), it yields 

1 f el2l / C*A°°e 2 C*^ 00 ^ C*a°°£ 2 \ /•* r .„ _ l2l , 

2 L M + (* 2 2 2T" J /„ I 1 ^'1'^ 

<^13 + ^/ / (|^H 2 +. 2 |W 3 | 2 ) 

2 jo JQ £ 

(63) 

where C13 depends on the bounded terms of r.h.s of ( |62|) . Differentiating now 
(TTTT) with respect to time and then testing the result with dt{u\ — uf), we get 

1 f \d t ul\ 2 dx + d 3 f f \Vd t u 3 £ \ 2 dxdr 

2 J Q £ Jo Jq e 

+ f I ( \(d t d 3 \Vu 3 £ \ 2 - (d t d t d 3 )\Vu 3 
Jo Jn? \2 



/ \Vu 3 £ \ 2 dxdT + ^- / \Vd t u 3 £ \ 2 dxdr 
2 Jo Jn? 2 Jo Jn e 

+ df + D 3 ft f ,„ n n „, , . eA 



[ [ \Vd t u 3 D \ 2 dxdT + - / / d t \u 3 £ \ 2 dxdr 

Jo Jn E 2 Jo Jr™ a 

nrA°° ft f 
\d t u 3 £ \ 2 dxdr + / / (|m 3 £ | 2 + \d t u 3 D \ 2 )dxdr 
jwa 2 jo Jr™ a 

+ ^ t I {\d t u 3 £ \ 2 + \d t u 3 D \ 2 )dxdr 
2 Jo Jr™ a 

fB°° rt r 

+ -S- / / (K| 2 + \d t u 3 £ \ 2 + \u 2 £ \ 2 + \d t u 3 D \ 2 )dxdr 



2 jo Jr™ 1 
Fh°° rt r 1 

+ -T- / / (tI^I 2 + 5|<W| 2 + |cW| 2 + |W|>ir. 
2 jo Jri" a 



Using (H) to deal with the boundary terms, we obtain 



1 r ( d°° C 1 *f'^ \ ft f 

- \d t u £ 3 \ 2 dx+(d 3 -^ — (3a°° + B°° + b°° + a°°5) / / \Vd t u 3 £ \ 2 dxdr 

2 J r2 £ y z z y Jo Jn s 
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<C U + C 15 [ [ \d t u 3 £ \ 2 dxdr (64) 
Jo Ju s 

+ C*b°° f f {\d t u 3 £ \ 2 + e 2 \Vd t u 3 £ \ 2 )dxdr (65) 
Jo Jn^ 

Adding ( 1631) and ( 164"1) and using ( 159]) and ( 160]) to get the desired result. 



4-1 Extension step 

Since we deal here with an oscillating system posed in a perforated domain, 
the natural next step is to extend all concentrations to the whole Q. We do 
this by following a two-steps procedure: In Step 1, we rely on the standard 
extension results indicated in section 14.21 to extend all active concentrations 
u\ (£ G {1, . . . ,4}) to Q. In step 2, we unfold the ode for u\ such that the 
unfolded concentration is defined on the fixed boundary T; see section 15.11 



4-2 Extension lemmas 

Since all the concentrations are defined in Q £ and f2f, to get macroscopic 
equations we need to extend them into Q. 

Remark 6 Take (f £ G L 2 (0,T; H 1 ^ 6 )) . Note that since our microscopic ge- 
ometry is sufficiently regular, we can speak in terms of extensions. Recall the 
linearity of the extension operator 

V £ : L 2 (0, T; H\Q £ )) -> L 2 {0, T; H\Q)) 

defined by V £ <p £ = (p £ . To keep notation simple, we denote the extension (p £ 
again by (p £ . 

Lemma 7 (Extension) Consider the geometry described in Section l2~l\ There 
exists an extension u £ of u £ such that 

(1) || u £ \\ LHY )< C || u £ \\ L2{YW) , foru £ G L 2 (Y W ) 

(2) || W \\ L *(Y)< C II W \\ L 2 {YW) , forWu £ G L 2 {Y W ) 

(3) || u £ || H i ( n)< C || u £ foru £ G 

Proof. For the proof of this Lemma, see Section 2 in [20J or compare Lemma 
5, p.214 in [30]. 

Definition 8 (Two-scale convergence cf. 131^32^ ) Let {u 6 } be a sequence of 
functions in L 2 ((0,T) x Q) (Q being an open set of R N ) where e being a 
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sequence of strictly positive numbers that tends to zero. {u £ } is said to two- 
scale converge to a unique function u (t, x, y) eI 2 ((0,T)x(lx Y) if and only 
if for any if) G C °°((0,T) x Q,Cg(Y)), we have 

lim £ ^o / / u £ ip(t,x, —)dxdt = / / uo(t,x,y)ip(t,x,y)dydxdt. (66) 
Jo Jn e Jn Jy 

We denote (6b}) by u £ u . 

Theorem 9 (i) From each bounded sequence {u £ } in L 2 ((0, T) x Q), one can 
extract a subsequence which two-scale converges to uo(t, x, y) G £ 2 ((0, T) x 
VI x Y). 

(ii) Let {u £ } be a bounded sequence in -ff 1 ((0, T) x Q), which converges weakly 
to a limit function u (t,x,y) G H l ((0, T) xflxY). Then there exists u G 
L 2 (f2; H^(Y)/M) such that up to a subsequence {u £ } two-scale converges 

to u (t, x, y) and Vu £ — V x u + V y u. 
(Hi) Let {u £ } and {eVu £ } be bounded sequences in L 2 ((0,T) x Q), then there 
existsuo G L 2 ((0,T) x Q; H^(Y)) such that up to a subsequence u £ and 
eVu £ two-scale converge to uo(t,x,y) and V y uo(t,x,y) respectively. 

Definition 10 (Two-scale convergence for e— periodic hypersurfaces fSjJlj) A 
sequence of functions {u £ } in L 2 ((0,T) x T £ ) is said to two-scale converge to a 
limit u G L 2 ((0, Tjxflxr) if and only if for any ip G C °°((0, T) x fi, C£(T)) 
we have 

lim £ ^Q£ [ [ u £ ip(t, x,—)da x dt = [ [ u (t,x,y)i/j(t,x,y)da y dxdt. 
Jo Jr e e Jn Jr 

Theorem 11 (i) From each bounded sequence {u 6 } G L 2 ((0,T) x r e ), one 
can extract a subsequence u e which two-scale converges to a function u G 
L 2 ((0,T) xflxT). 

(ii) If a sequence of functions {u £ } is bounded in L°°((0,T) x r e ), then u £ 
two-scale converges to a function uq G L°°((0, T) x Q x T). 

Proof. For proof of (i), see [33] and the one for (ii), see [To] . 

Lemma 12 Assume the hypotheses of Lemma [3] and Lemma [?] to hold. The 
a priori estimates lead to the following convergence results: 

(a) u £ -± Ui in L 2 (0, T; H^tt) for all i € {1,2,3, 4}, 

(b) u £ ^Ui inL°°({0,T) xty, 

(c) d t u £ -± d t Ui in L 2 ((0,T) x Q), 

(d) u £ ->■ ui in L 2 (0, T; H^iVL)) for \ < /3 < 1, also || u\ - u { ||L2(( 0iT)xre )->- 
as e — > 0, 

(e) u\ A u u Vu £ A V xUi + V y u a , u a G L 2 ((0, T) x Q; H^(Y)/R), 

(f) u\ A u 5; and u 5 G L°°((0,T) x x T sw ), 

(g) d t u% A d t u 5 , and u 5 G L 2 ((0,T) x x T sw ). 
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Proof, (a) and (b) are obtained as a direct consequence of the fact that u\ 
is bounded in L 2 (0, T; H 1 ^)) n L°°((0,T) x ft); up to a subsequence (still 
denoted by uf) u\ converges weakly to u { in L 2 (0,T; H 1 ^)) nL°°((0,T) x J]). 
A similar argument gives (c). To get (d), we use the compact embedding 
H p \Q) ^ H fi (n), for G (§, 1) and < < 0' < 1 (since ft has Lipschitz 
boundary). We have 

W := {ui G L 2 (0, T; H 1 ^)) and G L 2 ((0, T) x ft) for all i G {1, 2, 3, 4}} 

For a fixed e, W is compactly embedded in L 2 (0, T; H^(Q)) by the Lions- Aubin 
Lemma; cf. e.g. [34J. Using the trace inequality (jSJ) 

II «i - ||i2((o,T)xr 6 ) < Co II M f _ M « IU 2 (0,T;H^(n s ))) 
< C II M i - M « IU 2 (0,T;H^(n))) 

where || u\ — U; L \\l 2 (o,T;H^(Q))~ * as e -> 0. To investigate (e), (f) and (g), 
we use the notion of two-scale convergence as indicated in Definition [5] and 
[TU1 Since u\ are bounded in L 2 (0, T; if 1 (ft), up to a subsequence wf Mj in 
L 2 ((0, T) x ft x y), and Vitf A + V^, G L 2 ((0, T) x ft; #J(y)/R). 
By Theorem [HI w| in L°°((0, TjxOxT) converges two-scale to U5 in the same 
space and d t u% converges two-scale to d t u 5 in L 2 ((0,T) x ft x T). Due to the 
presence of the non-linear reaction rate on the interface T s e w , the convergences 
listed in Lemma [12] are still not sufficient to pass to the limit e — > in the 
microscopic model. To be more precise, we can pass to e — > in the pde's, 
but not in the ode. 



4-3 Cell problems 



In order to be able to formulate the upscaled equations, we define two classes 
of cell problems very much in the spirit of [9]. One class of problems will refer 
to the water-filled parts of the pore, while the second class will refer to the 
air-filled part of the pores. 

Definition 13 (Cell problems) The cell problems in water-filled part are given 
by 

3 

-V y .(D e (t,y)V yX i) = 1 £ d y k D iki(-t,y), ™ Y w , 

k=l 

By- JL 

-D £ (t, y)-^ = H D m (t, y)n k on T wa , 

k=l 
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for all i, £ G {1,2, 4} and Xi are Y-periodic in y. The cell problems in air-filled part 
are given by 

3 

-V y .(D 3 {t,y)V y ^) = Y,d yk D 3ki {t,y), m Y\ 

k=l 

-D 3 (t } y)^- = j2D 3ki (t,y)n k onT wa , 

k=l 

-D 3 (t, y)^ = jl D 3ki (t, V)n k on dY a - T wa , 

k=l 

for all i G {1,2,3} and q are Y-periodic in y. 



5 Two-scale limit equations 



Theorem 14 The sequences of the solutions of the weak formulation (fPj)-( fTgl) 

converges to the weak solution Ui,i G {1,2,3,4,5} as e — > such that Ui G 
H\0, T; L 2 (n)) n L 2 (0, T; H\tt)) n L°°((0, T) x Q) and u 5 G H\0, T; L 2 {Vt x 
r))nL°°((0, T) x OxT)). The weak formulation of the two-scale limit equations 
is given by 

ndtUi(t,x)(J)i(t,x)dxdt+ / / di(t)Vui(t, x)V(f)idxdt (67) 
_2 jo ./n 



/ / Fi(u)(j)idxdt for alii G {1,2,3,4}, 
Jo Jn 



where 



Fi(u) := -ki(t)ui(t, x) + k 2 (t)u 2 (t, x) 

r k 3 (t, x))Q(u 5 (t, x, y))da y , 



_ |Y|ir 

F 2 (m) := ki(t)ui(t, x) - k 2 (t)u 2 (t, x) + a(t)u 3 (t, x) - b(t)u 2 (t, x) 
~a(t)u 3 (t, x) + b(t)u 2 (t, x), 



F 4 



wzt/i i/ie initial values Ui(0,x) = Ui (x) for x <E Q, and 

rT 

d t u 5 (t, x, y)(fi 5 (t, x, y)dtdxda y 



y 

T ( ( 

= / / k 3 (t,y)R(u 1 (t,x))Q(u 5 (t,x,y))(j) 5 (t,x,y)dtdxda y , 
Jo Jnxr 

with u 5 (0,x,y) = u 5o (x,y) for x G Q, y G T sw . Also : = (0i, 2 , 03, 00 
[C°°({0,T) x fi)] 4 ; if, := tyuMaM e [C°°((0,T) x fi);C^(y)] 4 ; 
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k j {t)-=^rJ Y k j {t,y)dy, j G {1,2}, (69) 
^) := 7^7t/ a(t,y)dcry, (70) 



\Y\Jr 

b(t):=-j-^ b(t,y)d<j y , (71) 



3 r 

d Bj := Y, ( d «j(*> v) + duk(t, y)(8 ne d Vk Xj + 5 3e d yk <; j )dy, (72) 
k=i JY 

£e {1,2,3,4}, nG {1,2,4} 

with Xji^j being solutions of the cell problems defined in Definition [73, while 
S denotes here the Kronecker's symbol. 



Proof. We apply two-scale convergence techniques together with Lemma [T2l to 
get macroscopic equations. We take test functions incorporating the following 
oscillating behavior tpfo, x) = <fii(t, x) + eipi(t, x, f ), 0* G C°°((0, T) x fi), 0j e 
(7°°((0, T) xO, ; C# (K)), z G {1,2, 3, 4}. Applying^wo-scale convergence yields 



1*1 / / d t u i (p i (t,x)dxdt+ [[[ di(t,y)(V x Ui(t,x) 
Jo Jn Jo JqJy 

x 

+ V y Ui(t,x,y))(y x (pi(t,x) + V y ipi(t,x, -))dydxdt 

= / / fi(u)(f>i(t,x)dxdt. (73) 
jo Jn 



nfi(u) <pi(t, x)dxdt = — lim / / k\u\ (4>i(t, x) + eipUt, x, — ))dxdt 
A e^OJo Jn? e 



rT 



+ lim / / k^ui^ipiit, x) + sipAt, x, — ))dxdt 
e^o Jo J e 

— lime / rj(R(u\), Q(ul))(4>i(t, x) + eipUt, x, — ))da x dt. 
Jo Jrfu e 

Using Lemma Q2i we have 

f 1 (u)(j) 1 (t,x)dxdt = - / / / ki(t,y)ui(t, x)4>i(t, x)dydxdt 



/o ./<> jo Jnjy 

+ / / / k 2 {t 1 y)u2{t 1 x)4>i(t, x)dydxdt 
Jo Jn Jy 



nx 
d t u £ J(j)i(t, x) + stpi(t, x, —))da x dt. 
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nfi{u)4>i{t 1 x)dxdt = —\Y\ I 1 ki(t)ui(t, x)4>i(t, x)dxdt 
a Jo Jn 

+ 1*1 / / k 2 (t)u 2 (t,x)(j)i(t,x)dxdt 
Jo Jn 

— / / d t ur,(f)At,x)da v dxdt. 
Jo Jn Jr 



nf2(u)4>2(t,x)dxdt = \im / k\u\ (fait, x) + s^ 2 (t, x, — ))dxdt 

A e^OJo Jne £ 

— lim / / knuUfait, x) + efait, x, — ))dxdt 

e^OJo Jq,e £ 



T r X, 



+ lime / / a £ uU(j)2(t,x) + eip2(t,x, —))da x dt 

e->o Jo Jr™ a e 

— lime / / beUoifaityX) + eip2(t,x, —))da x dt. 

e->o Jo Jv™ a e 



nf 2 (u)(p2(t,x)dxdt—\Y\ [ [ k 1 (t)ui(t,x)4> 2 (t,x)dxdt 

a Jo Jn 

— \Y\ / k 2 {t)u2{t 1 x)(f) 2 (t, x)dxdt 

Jo Jn 

+ \Y\ / / a(t)us(t, x)(f) 2 (t, x)dxdt 

Jo Jn 

— \Y\ / b{t)u 2 {t 1 x)4> 2 (t, x)dxdt. 

Jo Jn 

We also have 

nf 3 (u)(f) 3 (t,x)dxdt = —\Y\ / / d(t)u 3 (t, x)(f) 3 (t, x)dxdt 
a Jo Jn 

+ \Y\ I [ b(t)u 2 (t,x)(f)3(t,x)dxdt 

Jo Jn 

and 



nfi{u)(t) i {t 1 x)dxdt=\Y\ I I k\{t)ui{t, x)0 4 (t, x)dxdt. 
a Jo Jn 

We set 4>i = 0, i G {1, 2, 3, 4} in ( 1731) to calculate the expression of the known 
function u\ and obtain 

m di(t,y)(\7 x iii(t,x) + S7yUi{t,x,y))S7 y ipi{t } x } -)dydxdt = 0, forall ^. 
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Since u\ depends linearly on V x Ui, it can be denned as 

3 

Ui := ^2d Xj Ui(5 in Xj{t,y) + S^j(t,y)) for n G {1,2,4} 

where the function Xji^j are the unique solutions of the cell problems defined 
in Definition ITBI Setting ipi = in ([75]) . we get 



^ f f it d ijk(t,y)( d x k Ui(t, 



3 

+ X! ( 5 in^ fc Xm + hid Vh S m )d Xm Ui{t,x))d Xj (t> k {t,x)dydxdt 

m=l 

= 1^1 / / dijkd Xk Ui(t, x)d x J>i(t, x)dxdt. 
Jo Jn j k=1 

Hence, the coefficients (entering the effective diffusion tensor) are given by 
1 3 r 

dijk ■= TTTT / ( rf «i(*> + d «*(*> y)( s indy k Xj + hidy k ^j)dy. (74) 
l r I fc=l ,/y 

j G {1,2,3, 4}, n G {1, 2, 4} and j, fc G {1, 2, 3}. 



5. i Passing to the limit e — > m (T7 



It is not yet possible to pass to the limit e — > with the convergence results 
stated in Lemma [T2l To overcome this difficulty, we use the notion of periodic 
unfolding. It si worth mentioning that there is an intimate link between the 
two-scale convergence and weak convergence of the unfolded sequences; see 
[35p5] . The key idea is: Instead of getting strong convergence for it|, obtain 
strong convergence for the periodic unfolding of u\. 

Definition 15 For e > 0, the boundary unfolding of a measurable function ip 
posed on oscillating surface T £ is defined by 

T*ip(x,y) = ip(ey + ek), yeT,xett 

where k := [|] denotes the unique integer combination H^ =1 kjej of the periods 
such that x— [-] belongs to Y . Note that the oscillation due to the perforations 
are shifted into the second variable y which belongs to fixed surface V . 

Lemma 16 If u £ converges two-scale to u and T§u e converges weakly to u* 
in L 2 ((0,T) x fi;L|(r)), then u = u* a.e. m (0,T) xllxf. 
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Proof. The proof details for this statement can be found in Lemma 4.6 of [T5] . 



Lemma 17 If (p £ L 2 ((0,T) x r e ) ; i/ien £/ie following identity holds 
1 
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K&VlU 2 ((o,T)xnxr) — e||y|U 2 ((o,T)xr e )- 



Proof. Consider ^ 

ivf ^l^nxr) = Tv\ L ^ \ T bf\ 2 dxda y = — T^ 2 dxda y , 

= i^ s Li / / Tfa 2 dxdo- y = i^t^Li /, d.*c i ^ 2 da y , 
| r | J £ (y+fc) | Y | J £ (y+A:) Jr 

= Yl =1 e* J^ 2 da y . 
Changing variable z = e(y + k), where k = [-], we get 



TTFj\ T bV\h (nxr) = s Li^ 3 / <^ 2 ^y = s Li^ / , = e / ¥> 2 <*<t z . 

|r| v y Jr J e (r+fc) Jr^ 

This completes the proof of (1170 . 



Lemma 18 IftpE L 2 {VL), then T^(p —>■</? as e — >• strongly in L 2 {Vt x T). 
Proof. See in [36f3T] for proof details. 



Using the boundary unfolding operator T 6 e , we unfold the ode ffT3]) . Changing 
the variable, x = ey + ek (for x £ r| w ) to the fixed domain (0, T) x fi x T, we 
have 



d t T^ul(t,x,y) = r,(T^ul(t,x,y),T^ul(t,x,y)). (75) 

In the remainder of this section, we prove that T£u £ 5 converges strongly to in 
L 2 (Q x T). From the two-scale convergence of w|, we obtain weak convergence 
of T e u e 5 to u 5 in L°°((0, T) x Q; Z^ er (r)). We start with showing that {T b e w!} is 
a Cauchy sequence in L 2 (Q xT). To this end, we choose m, n £ N with n > m 
arbitrary. Writing down (!75l) for the two different choices of e (i.e. = e n and 
£j — £ m), we obtain after subtracting the corresponding equations that 
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= / - klR{T^ur)Q{T^ur)) 

Jnxr 

x (T^ l u e 5 n - T^ m u £ 5 m )da y dx, 

0°° r 
< fefcaC^- + CQSwpaxiK'X"!) / |T 6 e »< - T^ut m \ 2 da y dx 

Z JflxF 

+ 3 * V / |2£"tt?» - T^ttM 2 ^^. (76) 

To get (175"|) . we have used the uniform boundedness of T^ n u\ n . We consider 
now 



/ m n u?--Tf n ur\ 2 da v dx 
Jnxr 

< J v {\n n u? - + m n ui - Ul \ 2 )da y dx 

+ / (\T* m Ul - u,\ 2 + \T™u\ m - T™ Ul \ 2 )da y dx. (77) 
Jnxr 

Since u\ is constant w.r.t. y, we have that T^ m u\ — > u\ strongly in L 2 ((0,T) x 
Q x T) as e — >■ 0. From Lemma fT7l we conclude that 

/ I^Wi — T^ui\ 2 da y dx < e \ \u\ — ui\ 2 da y dx < eC. 
Jnxr Jr s 

(1771) turns out to be 



/ | T^u\ n - T^ n u\ m \ 2 da v dxdt <C(e n + e m ), 
Jnxr 

while (17611 becomes 

d t [ \T™ul n - T™ul m \ 2 da y dx < Cis / \Ttu e 5 n - T^ul m \ 2 da y dx + — , 

where C 15 := kfc R {^ + c Q su P nxv\T^u\ n \) and C 16 := ^f^C. The Gron- 
wall's inequality gives 

II rpen n ,en rpem^em \\ ^ ^16 /'VoN 

1Mb % M 5 IU 2 (Qxr)< — • (78) 

By (175|) . {T^Mg} is a Cauchy sequence. Now, we take the two-scale limit in the 
ode ( T75|) to get 



nx f x 

d t T^ul(j)i(t, x, — )da x dt = lime / / ?7(T fe e w?, T fe e -u?)0i(t, x, —)da x dt. 



Consequently, we have 
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rT r 

d t u 5 (p 5 (t, x, y)dxda y dt 



rT 



= lime / / T^R(T*ul)Q(u £ 5 )<p 5 {t,x,-)da x dt, 

-lime / T / T?kiR(T?ul)Q(u 5 U 5 (t,x,-)dcr x dt 

+ lime / T / T^lR{T^u\){Q{T^ul) - Q(u 5 ))(f> 5 (t,x, -)da x dt. 
e->0 Jo Jy% w e 

(79) 

By (A2) and the strong convergence of u\, the first term on the right hand 
side of (1791 converges two- scale to 

III k 3 (t,y)R(u 1 )Q(u 5 )(j) 5 (t,x,y)da y dxdt, 
Jo Ju Jr sw 

while the second integral of (179"]) 

e f T I T^klR{T^u\){Q{T^ul) - Q(u 5 ))(/> 5 (t,x, -)da x dt 
Jo Jv% w e 

i 

< e ( f T [ \TZklR(T b e ul)ci> 5 (t,x, -)\ 2 da x dt) ' ■ 
\Jo Jri w e / 



\Q(T^ut)-Q(u 5 )\ 2 da x dt)\ 

JVf w I 



->0 as e — >■ 0. 

At this point, we have used again (A2) in combination with the strong con- 
vergence of T^u\. So, as result of passing to the limit e — > in (p~3|) we get 

QSSD- 



It is worth noting that the weak solution to the two-scale model inherits a.e. 
the positivity and boundedness properties from the corresponding properties 
of the weak solution of the microscopic model. Now, it only remains to ensure 
the uniqueness of weak solutions to the upscaled model. 



Lemma 19 (Uniqueness of solutions of (fffy- ffity) Assume (Al)-(A6). There 
exists at most one weak solution to the two-scale limit problem (fFfy and fiUfy) . 

Proof. Suppose there are two weak solutions to the two-scale limit problem 
(ill, u 2i u i, u ly u i) with j e {1, 2}. We denote ue = u\ — uj, £ G {1, 2, 3, 4} and 
choose as test function (pi = ui. After straightforward calculations, we have 
from (J6ED 



rt 

\ul — ul\ <C Itii — u\\dr. 
Jo 
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Take <pi = u\ in (157)1 to obtain 



[ [ d t \ui\ 2 dxdt + di [ [ \Vui\ 2 dxdt 
Jo Jn Jo Jn 



2 Jo Jn Jo Jn 

<-h f I d t \u 1 \ 2 dxdt + ^- f ( Oi| 2 + \u 2 \ 2 )dxdt 
Jo Jn 2 Jo Jn 

k°° f t f 

+ t^tCrCqMx / / (u\- ul)u x dxda y dt 
\Y | " Jo Jnxr sw 

k°° ft c 

+ t^c r Q°° / \ Ul \ 2 dxda y dt. (81) 
\Y | Jo Jnxr s » 

Using ( )80|) together with the trace inequality for fixed domains, see section 
5.5 Theorem 1 in [38] and also the fact that U\ is independent of y in ( IHTj) . we 

get 



/ / d t \u 1 \ 2 dxdt + {2d 1 -k^c R C*{5M 1 + Q co )) f f {Vu^dxdr 
Jo Jn ' Jo Jn 

+ 2k~! f f dtlu^dxdr 
Jo Jn 

< (fc 2 °° + k^c R C*(6M 1 + Q°°)) f [ (\ Ul \ 2 + \u 2 \ 2 )dxdr 

Jo Jn 



io Jn 

+ k™5c R c Q M 1 C* [ [ [ {\uA 2 + \Vu 1 \ 2 )dsdxdT. 
Jo Jn Jo 

2di-fe|° c R C* 

V 3 



For suitable choice of 5 e]0, 1 k oo l R c*M 1 — [' we have 

/ / d t \ui\ 2 dxdt + d 1 [ [ \Vu 1 \ 2 dxdr + 2k 1 [ [ d t \ui\ 2 dxdT 
Jo Jn Jo Jn Jo Jn 



< (fc 2 °° + k^c R C*{5M x + Q°°)) [ T [ (K| 2 + \u 2 \ 2 )dxdr 

Jo Jn 



+ ^c R c Q M 1 C* ( [ f \\ Ul \ 2 + \V Ul \ 2 )dsdxdr. (82) 
o * Jo Jn Jo 

Take (fi 2 = u 2 in fl67j) . we get 

- / [ d t \u 2 \ 2 dxdr + d 2 [ [ \ Vu 2 \ 2 dxdr 
2 Jo Jn Jo Jn 



<-~k 2 f I d t \u 2 \ 2 dxdT + ^- f ( {\ Ul \ 2 + \u 2 \ 2 )dxdr 
Jo Jn 2 Jo Jn 

+ d°° / / u 2 u^dxdt — b / / \u 2 \ 2 dxdr. 
Jo Jn Jo Jn 



I I d t \u 2 \ 2 dxdr + d 2 f f \Vu 2 \ 2 dxdr 
Jo Jn Jo Jn 



< (k°° + a°°) f f (\ Ul \ 2 + \u 2 \ 2 + \u 3 \ 2 )dxdr. (83) 
Jo Jn 
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Similarly, we obtain from (|67|) 



f j d t \u 3 \ 2 dxdr + d 3 [ [ \Vu 3 \ 2 dxdr < ~b°° [ j (\u 2 \ 2 + \u 3 \ 2 )dxd($A) 
h Jn Jo Jn Jo Jn 



/* / d t \u 4 \ 2 dxdT + d 4 f [ \Vu 4 \ 2 dxdr <k<? T [ (M 2 + |w 3 | 2 )dx485) 
Jo Jn Jo Jn Jo Jn 

Adding side by side (1521 - (1551) and applying Gronwall's inequality to the cor- 
responding result, we receive 

E.i x / \ Ui \ 2 dx + effi^i / / \Vui\ 2 dxdr + d [ f \ui\ 2 dxdr<0. (86) 
Jn Jo Jn Jo Jn 

In flHBI) . we have d := min{di,d 2 ,d 3 ,d 4: ,ki} > 0. Taking in flSTl) supremum 
over (0,T), we obtain 

£*=i / H 2 dx + cxTti / T / |Vxx;| 2 dxdr < 0, (87) 
in jo in 

which concludes the proof of the Lemma. 

Lemma 20 (Strong formulation of the two-scale limit equations) Assume the 
hypothesis of Lemma [HI to hold. Then the strong formulation of the two-scale 
limit equations (for all t G (0, T) ) reads 

d t ui(t, x) + V • (— diVui(t, x)) 

= —ki(t)ui(t, x) + k 2 (t)u 2 (t, x) 

— 7777 / k 3 (t,y)R{u 1 {t,x))Q{u 5 {t,x,y))da y , x G Q 
y Jr aw 

Ui(0,x) = u w (x), x G 0, 
n ■ (— diVixi(i, x)) =0, x G <9f2 

(88) 

d t u 2 (t, x) + V ■ {-d 2 Vu 2 (t, x)) = ki(t)ux(t, x) - k 2 (t)u 2 (t, x) 

+ d(t)u 3 (t, x) - b(t)u 2 (t, x), x efl, , . 

ix 2 (0,x) = li20(^), z G i2, 
n ■ (—d 2 Vu 2 (t,x)) = 0, x G <9f2, 

d t u 3 (t, x) + V • (~d 3 Vu 3 (t, x)) = —a(t)u s (t, x) + b(t)u 2 (t, x), x G fi, 
ix 3 (0,x) = u 30 (x), x G fi, 

1x3 (t, x) = (x), x g r D , 

rx ■ (-d 3 Vxx 3 (t, x)) = 0, x G 

(90) 
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dtUi(t, x) + V • (— diVuiit, x)) = ki(t)ui(t, x), x G Q, 
1*4(0, x) = U4o(x), (91) 
n ■ (—d^Vu^tjx)) = 0, x G <9f2 

<9 4 u 5 (t, x, 2/) = k 3 (t, ?/)i?(M 1 (t 1 a;))Q(u5(t, x, y)), x ett,y E T sw , 
m 5 (0, x, y) = u 50 {x, y) x ett,y e T sw , 



where di,i G {1,2,3,4} and kj,j G {1,2} are defined in Theorem 14 
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